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, In this article, we formulate the study of the unitary time evolution of systems consisting 

of an infinite number of uncoupled time-dependent harmonic oscillators in mathematically 
rigorous terms. We base this analysis on the theory of a single one-dimensional time- 
dependent oscillator, for which we first summarize some basic results concerning the unitary 
, implementability of the dynamics. This is done by employing techniques different from those 

used so far to derive the Feynman propagator. In particular, we calculate the transition 
(— I , amplitudes for the usual harmonic oscillator eigenstates and define suitable semiclassical 

states for some physically relevant models. We then explore the possible extension of this 
study to infinite dimensional dynamical systems. Specifically, we construct Schrodinger 
functional representations in terms of appropriate probability spaces, analyze the unitarity 
of the time evolution, and probe the existence of semiclassical states for a wide range of 

■ physical systems, particularly, the well-known Minkowskian free scalar fields and Gowdy 

^ \ cosmological models. 
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I. INTRODUCTION 

The quantum time-dependent harmonic oscillator (TDHO) has received a lot of attention due 
to its usefulness to describe the dynamics of many physical systems. This is the case, for example, 
of radiation fields propagating outside time-dependent laser sources or in spatial regions filled with 
time-dependent dielectric constant matter [ij. The behavior of ions in Paul traps can also 

be described by one-dimensional harmonic oscillators with time-dependent frequencies. The math- 
ematical aspects of the quantum TDHO and its applications to more general theoretical models 
have been profusely analyzed in the literature [1, [aQjS]. In particular, they have been studied in 
the context of the search of exact invariants for nonstationary quantum systems. This method was 
introduced for the first time by Lewis and Riesenfeld and proved to be especially useful to 

generate exact solutions to the Schrodinger equation and also to probe the existence and properties 
of semiclassical states for these systems. 
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In the context of infinite dimensional dynamical systems, physical theoretical models with in- 
finitely many time-dependent oscillators naturally appear in quantum field theory in curved space- 
times 111, [l^ and also in the reduced phase space description of some midisuperspace models in 
general relativity 13, 14, iB, 3]. It is clear that, in the quest for a suitable quantization for these 
systems, the understanding of the special features of the single TDHO is particularly advisable. 
This fact motivates us to summarize in a rigorous and self-contained way the theory of the quan- 
tum TDHO in the first part of the paper by (implicitly) making use of the theory of invariants. 
Although some of the basic results such as the formula of the Feynman propagator are certainly 
well-known, they will be recovered by using novel techniques, and contrasted with the expressions 
obtained in the existing literature. This procedure will facilitate the generalization to field theories 
in the second part of the paper, where we apply the ideas exposed for the one-dimensional oscillator 
in order to obtain the propagator for infinite-dimensional systems and discuss applications both to 
quantum field theory and quantum gravity. The structure of the article is the following. 

In section\l^ we analyze some relevant properties of the classical TDHO equation, in particu- 
lar, its connection with the so-called Ermakov-Pinney (EP) equation 171. llSl. Il9l|. which plays an 
auxiliary role in the calculation of invariants for nonquadratic Hamiltonian systems Many of 
the quantities derived in this section are used afterward to obtain a simple and closed expression 
for the unitary evolution operator of the quantum TDHO. 

In sectioning we first briefiy review the definition of the abstract Weyl C*-algebra of quantum 
observables for the TDHO, the uniqueness of all regular irreducible representations of the canonical 
commutation relations, and the unitary implementability of the symplectic transformations that 
characterize the classical time evolution. Once a concrete representation is fixed, we construct the 
unitary evolution operator by introducing some suitable displacement and squeeze operators j2l| . 
In our discussion, the use of the auxiliary EP equation is appropriately interpreted as a natural 
way to solve the Schrodinger equation and obtain an expression of the evolution operator valid for 
all values of the time parameter. We then show that this method is especially useful to derive the 
Feynman propagator, whose calculation follows readily in this context. We obtain an expression in 
agreement with those previously derived in the existing literature, where (more complicated) path- 
integration techniques are often employed [i^, ■ We also analyze in this section the calculation 
of the transition amplitudes for the usual harmonic oscillator eigenstates and, as a particular case, 
the instability of the vacuum state as a direct consequence of the nonautonomous nature of the 
system. 

The eigenstates of the Lewis invariant provide a family of state vectors closed under time 
evolution, depending on a particular solution to the EP equation, that generalize the minimal wave 
packets of the harmonic oscillator with constant frequency. These states, however, do not have the 
usual properties of the ordinary coherent states -not even the ones associated with the squeezed 
states- and can be taken as semiclassical states just in case they are well-behaved enough. In sec- 
tion \IV[ we analyze the construction of semiclassical states for some physically relevant systems, 
such as the vertically driven pendulum and, particularly, the class of TDHO equations that occur in 
the well-known Gowdy cosmological models 2j] , as a previous step to generalize this construction 
to field theories. 

SectionlM is precisely devoted to the extension of the previous study to linear dynamical sys- 
tems with infinite degrees of freedom governed by nonautonomous quantum Hamiltonians that 
can be interpreted as systems of infinite uncoupled harmonic oscillators with time-dependent fre- 
quencies. This fact allows us to give a straightforward procedure to obtain the unitary evolution 
operator, following the discussion developed for a single oscillator. We particularize our results 
to the well-known Minkowskian free scalar fields and also to the Gowdy cosmologies, that have 
attracted considerable attention in recent years as appealing frameworks to test quantum gravity 
theories (see 14, 15, 1^, 25, 26, 27 1 and references therein). Making use of Schrodinger representa- 
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tions, where states act as functionals on appropriate quantum configuration spaces, we construct 
the analog of the one-dimensional propagator. We also discuss the difficulties that arise when 
dealing with infinite dimensional systems -specifically, the impossibility of unitarily implementing 
some symplectic transformations- and their implications for the search of semiclassical states. We 
conclude the paper with some final comments and remarks in section and appendix HI 

Throughout the paper, we will take units such that the Planck constant h, the light velocity c, 
and the characteristic mass of the system under study are equal to one. For any z £ C \ (—00, 0] , 
^/z will denote the unique square root of z such that Ke{^/z) is strictly positive. 



II. PROPERTIES OF THE TDHO EQUATION 

We will review in this section some properties of the classical equation of motion of a single 
harmonic oscillator with time-dependent frequency, from now on referred to as the TDHO equation, 
given by 

u{t) + K{t)u{t) = , t e / = C M, (2.1) 

where k : / — > M is a real- valued continuous function and time-derivatives are denoted by dots. 
Given an initial time £ I , let Ct^ and s^g be the independent solutions of (j2.1|) such that 
cto(^o) = stoC^o) = 1 and stp(to) = ctg(to) = 0. These can be written in terms of any set of 
independent solutions to (|2.ip . say ui and U2, as 

U2(to)ui{t) - Ui{to)u2{t) Ui{to)u2{t) - U2{to)ui{t) 

Cto [t) = TTTf ^ > ■Sio [t) = TTTf ^ , 2.2) 

W[Ui,U2) W[Ui,U2) 

where {to,t) £ I x I and W{ui,U2) := ^1112 — U1U2 denotes the (time-independent) Wronskian of 
ui and U2- In what follows, we will use the notation c{t,to) := Cfo(t), c(t,to) '■= cto{t), s{t,to) := 
stQ{t), and s(t, ^o) '■= ^toii)- Note that the s function belongs to the class C^(I x /) , whereas 
c{-,to) G C^il) and c(t, •) E C^(I) . As a concrete example, for the time-independent harmonic 
oscillator (TIHO) with constant frequency K{t) = kq £M, we simply get (w > 0) 

Ko = , c{t,to) = cos{{t — tQ)uj) , s{t,tQ) = uj^^ sin{{t — tQ)uj) ; (2-3) 
Ko = 0, c{t,to) = l, s{t,to) = t-to; (2.4) 

Ko = — , c(t, to) = cosh((t — to)"^) 1 s{t, to) = uj^^ sinh{{t — tQ)uj) . (2-5) 

In fact, as well known from Sturm's theory, the c and s functions corresponding to arbitrary fre- 
quencies share several properties with the usual cosine and sine functions. Firstly, their Wronskian 
is normalized to unit, W{c, s) = 1. Hence, if one of them vanishes for some time t = t*, then the 
other is automatically different from zero at that instant. In view of this condition and Eq. ()2.2p . 
their time-derivatives satisfy 

s{t, to) = c{to,t) , c(t, to) = — — , (2.6) 

s{t,to) 

where the last equation must be understood as a limit for those values of the time parameter 
t=K such that s(t=K,to) = 0. The odd character of the sine function translates into the condition 
s{to, t) = —s{t, to)- Finally, the well-known formula for the sine of a sum of angles can be generalized 
to 



S{t2,tl) = c{ti,to)s{t2,to) - c{t2,to)s{ti,to) 



(2.7) 
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It is well known that solutions to the TDHO equation (j2.ip are related to certain non-linear 
differential equations. Here, we will restrict our attention to the so-called Ermakov-Pinney (EP) 
equation (see 17, the interested reader is strongly advised to consult the historical account of 
[191 1 ^'^d references therein). Let 



A 



On ai2 

0.12 022 



be a positive definite quadratic form with det A = 1. Then, the (never vanishing) function p : / 
(0, +oo) defined as 



p{t) := ^/anc^^{t, to) + a22S^{t, to) + 2ai2s(t, to)c(t, to) (2.8) 
satisfies the EP equation 

p(t) + K(t)p(t) = , tel. (2.9) 

According to Eq. (j2.2p . the most general analytic solution to Eq. (j2.9p can be written as 2^, 29] 

pit) = \/bnuj{t) + b22ul{t) + 2buui{t)u2{t) , (2.10) 

where, as a consequence of (12. Sh and ()2.9p . the coefficients 612, 622 G 1^ satisfy Ty^(ni,U2) = 
(611622 — ^12)"^ > 0- Conversely, given any solution to the EP equation it is possible to find 
the general solution to the TDHO equation. Indeed, it is straightforward to prove the following 
theorem. 



Theorem 1. Let p be any solution to the EP equation i2.9\) : then, the c and s solutions to 12. 1\) 
are given by 

sit,to) = p(t)p(to) sin ^ , {t,to)eIxI. (2.12) 

Remark 1. By using Eq. (12. lip and ()2.12p . it is possible to find other p-independent objects. For 
example, the combination 

cos ( / -^tt) + P{to)p{t) sin ( / -Tr\] = "^(^o, i) = s{t, to) 



Pit) \kp'ir)J \JtoPHr) 
and the zeros of s(t,to)) characterized by 

dr 



to 



P'ir) 



(modvr) 



are independent of the particular solution p to the EP equation. These results will be profusely 
applied along the article. 
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III. UNITARY QUANTUM TIME EVOLUTION 
A. General framework 

The canonical phase space description of the classical system under consideration consists of 
a nonautonomous Hamiltonian system (/ x T,dt,u, H(t)). Here, T := denotes the space of 
Cauchy data {q,p) endowed with the usual symplectic structure i^{{qi,Pi), {q2,P2)) '■= P1Q2 — 
P2Qi, V (^2,^2) £ r. The triplet (/ x r,dt,a;) then has the mathematical structure of a 

cosymplectic vector space (see [s^] for more details). The time-dependent Hamiltonian H{t) : T — s- 
M, t G /, is given by 

H{t,q,p):=^(p^ + K{t)q^y (3.1) 

The solution to the corresponding Hamilton equations with initial Cauchy data {q,p) at time to 
can be written down as 

qH{t,to)\_^ f q\ ^ _ f c{t,to) s{t,to)\ .„„x 

PH{t, to) J - ■ U y ' V cit, to) sit, to) J ■ ^^-^^ 

Note that the properties stated in section O about the c and s solutions to the TDHO equation 
(f2T|) guarantee that 'T^t,to) ^ SL{2,R) = SP{1,R) for aU {t,to) £ I x I , i.e., the classical time 
evolution is implemented by symplectic transformations. 

We now formulate the quantum theory of the TDHO by defining an appropriate abstract C*- 



algebra of quantum observables 31[. This algebraic approach, although being more complicated 
than the traditional canonical quantization for this system, will facilitate the study of field theories 
in subsequent sections. We first realize that, as a consequence of the linearity of F, the set of 
elementary classical observables can be identified with the R-vector space generated by linear 
functionals on T. Every pair A := (—6, a) G F has an associated functional Fx : T — > M such that, 
for all X = {q,p) £ F, Fx{X) := a;(A, X) = aq + bp . In order to quantize the system, we introduce 
the abstract Weyl C*-algebra of quantum observables on T, '^{T), generated by the unitary Weyl 
operators W{\) = ex.p{iFx), X £ T, satisfying 

WiXi)* = Wi-Xi), WiXi)W{X2) = exp{iu{Xi,X2)/2)WiXi + X2), VAi,A2Gr. (3.3) 



According to von Neumann's uniqueness theorem [3ll . |32| |. all regular irreducible representations 
vr : WiT) ^{J^) of the Weyl C*-algebra into separable Hilbert spaces {J^, {■ \ •)) are uni- 
tarily equivalent. Here, ^{J^f) denotes the collection of bounded linear operators on J^. A 
*-homomorphism tt : W{T) — > ,'^{Jif) is said to be a regular irreducible representation if it has {0} 
and as the only closed vr-invariant subspaces, and tt(W{0, a)) and Tr{W{—b, 0)) are strongly con- 
tinuous in the a and b parameters, respectively. A well-known solution is given by the Schrodinger 
representation vTs : #'(r) — > =^(L^(R, d(?)) into the Hilbert space L^(]R, dg) where, for all pure 
states ip G L'^{R,dq), 

{7rs{W{X)) ■ tp){q) := exp {-iab/2) exp{iaq)ip{q + b) , A = (-6, a) G T . 

Thanks to the regularity condition, the usual Heisenberg algebra can be recovered in a definite 
sense from the Weyl C*-algebra. The strong continuity of ■Ks{W{0,a)) and TTs{W{—b,0)) in the 
real variables a and b ensures, by virtue of Stone's theorem, the existence of (unbounded) self- 
adjoint generators Q and P with dense domains in L^(R, dg). In particular, the Schwartz space 
^(R) of smooth rapidly decreasing functions in R is a common invariant dense domain of essential 
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self-adjointness for Q and P, where the usual Heisenberg algebra is satisfied. For all ip S =5^(M), 
we have {Qil^){q) = qip{q) and {Pip){q) = —i'il)'{q)^ where ip' denotes the derivative of '0- 

Another possibility is to represent the canonical commutation relations (CCR) in the space 
L^(M, d;Ua) where, given some a G C \ {0}, ^q, denotes the Gaussian probability measure 

V27r|a| V 2|a|^y 

To each a there corresponds a family of unitary transformations V^(/3) : d(7) L^(M, d/i^) 

connecting the standard Hilbert space with the new one in the form 

^[q) = = (^/2^|a|)'/'exp(-i^gV(2«))^('?), (3.4) 

where the complex numbers (3 must satisfy a/3 — /3a = i . Note that the unitary transformations 
VaiP) map the 'vacuum' state 

Mq) = (\^|a|)"^/^exp (i^gV(2a)) G L\R,dq) 

into the unit function ^'o('z) = {Va{(3)'4'o){q) = 1 G -^^^(1^, d/ia) . In these cases, the position and 
momentum operators act on state vectors as 

{Q^){q)=q^{q) and {P^){q) = -i^' (q) + ^q^{q) , 

a 

where, with the aim of simplifying the notation, Q and P respectively denote the Schrodinger trans- 
formed operators Va{(3)QVa{(3)~^ and Va{P)PVa{P)~^ with common dense domain Va{(3)y{E.) C 
L2(R,d/i„) . 

Any regular irreducible representation vr : W{T) — > ^(r^) is stable under time evolution, i.e., 
there exists a (biparametric) family of unitary operators U (t, to) : ■, the so-called quantum 

time evolution operator, such that 

C/-nt,to)vr(l^(A))C/(t,to) = ^(T^Cr(t,t„)A)), VVF(A) G #-(?) , V(t,to)e/x/, (3.5) 

with T^^t^f^-^ defined in Eq. (j3.2p . These relations determine C/(t,to) univocally up to phase. It is 
important to notice at this point that, if the classical evolution has singularities at the boundary 
of the interval /, they also occur for the quantum dynamics, i.e., there is no resolution of classical 
singularities. On the other hand, we will check in section [2 that the unitary implementability of 
symplectic transformations like those corresponding to the classical time evolution is not directly 
guaranteed for infinite-dimensional systems. The Heisenberg equations for Q and P can be solved 
just by the same expressions involved in the classical solutions ()3.2|) . i.e.. 

With more generality, given any well-behaved (analytic) classical observable F : F ^ M for the 
TDHO, the time evolution of its quantum counterpart Fnit^tQ) := U^^{t,to)F{Q, P)U{t,to) in 
the Heisenberg picture is simply given by 



FHit,to) = F{QHit,to),PH{t,to)) = F{c{t,to)Q + sit,to)P,c{t,to)Q + sit,to)P) . (3.7) 
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Hence, the matrix elements {^2\U-'^{t2,ti)F{Q,P)U{t2,ti)'^i), ^'i,^'2 G can be computed 
without the exphcit knowledge of the unitary evolution operator. This is also the case for the 
probability transitions Prob(^'2,t2 | = 1(^*2 I U{t2,ti)'ifi)\'^ , as will be discussed in detail in 

subsection UIIIA The commutators of time-evolved observables can be also calculated without the 
concrete expression of U{t2,ti). For instance, from Eq. (j3.7p we easily obtain 

[QHitl,to),QH{t2,to)] = is{ti,t2)l , 

where we have used the relation (j2.7p stated in section O As expected, the commutator given 
above is proportional to the identity operator and independent of the choice of the initial time to ■ 
Note, in contrast with the transition probabilities, that the calculation of transition amplitudes 
of the type (^^2 I ^^(^21^1)^1) does require the explicit knowledge of (the phase of) the evolution 
operator. This is also the case for the (strong) derivatives of both U{-,to) and U(t, •) . 

The dynamics of the quantum TDHO is governed by an (unbounded) nonautonomous Hamil- 
tonian operator H{t) : M' — > M', t £ I, satisfying 

Uit,to) = -iHit)U{t,to). (3.8) 



Given the quadratic nature of the classical Hamiltonian (jS.ip . H(t) must coincide with the operator 
directly promoted from the classical function modulo a t-dependent real term proportional to the 
identity 1 that encodes the election of C/(t,to) satisfying Eq. (jS.Sp . For a concrete representation 
of the CCR, we will simply take 

H{t) := i (P2 + ^(t)Q2) (3 9) 

This choice fixes U (t, to) uniquely. The Hamiltonian (jS.Op is a self-adjoint operator with dense 
domain &H{t) -equal to C^(M) in the standard Schrodinger representation- for each value of the 
time parameter t . We will prove the following theorem in the next subsections. 

Theorem 2. The action of the unitary TDHO evolution operator U{t,to) corresponding to the 
Hamiltonian i3. 9\) on any state vector £ C L^(M, dg) in the traditional Schrodinger 

representation is given by 



{U{t,to)^p){q)= [ K{q,t-qo,to)^P{qo 

Jm. 



dqo, 



where the propagator K{q,t;qo,tQ) depends on the times to o,nd t through the classical TDHO 
solutions c and s . Explicitly, 

K{q, t; go, to) = s"'/'(i, to) exp (^^^^ {c{to,t)q^ + c{t, to)ql - 2qqo)^ , (3.10) 

wherever s{t,to) ^ 0, and 

Kiq,t;qo,to)=c-^/\t,to)exp (^^^) %0 - 'zMt, to)) (3.11) 

ifs{t,to)=0. 

Remark 2. Given a solution u{t) to the TDHO equation (|2.ip which is positive in some interval 
(to, to + e) C /, e > 0, we define 

u^{t,to) := exp (^ie7Txn{u;t,to))\u{t)\^ , eGM, tel, 
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where m(n;t,to) E Z is the index function of u, with m{u;to,to) = 0, defined in such a way that 
m(n; t2,to) — m{u] ti,to), ti < t2, gives the number of zeros of «(•, to) in the interval (ti, • Finally, 
6{q) denotes the Dirac delta distribution. 

Remark 3. Let -i? : I ^ R be a real-valued continuous function and consider the Hamiltonian 

Hi{t) := H{t)+^{t)l 



defined in terms of (|3.9p . The unitary evolution Ui{t, to) associated with Hi{t) satisfying Eq. (jS.Sp 
gives rise to the propagator 

Ki{q,t;qo,to) = K{q,t : qo,to)e^p (^-i i9(r)dr^ . (3.12) 

Note that t/f ^(t, to) O Ui{t, to) = U^'^it, to) O U{t, to) for any quantum observable O. 

Remark 4. In the -L^(M, d/XQ,)-representation defined by the unitary transformation Va{(3) (see 
Eq. (j3.4p ). the evolution is given by 

{U{t,to)^)iq) = / Ka/3(9,t;go,to)^(go)d/ia(go), 

JR 

where 

f i/3 ij3 \ 

Ka/3{q, t; go, to) := V2Tr\a\ exp f —q^ - —q^ j K{q, t; go, to) • (3.13) 

B. Constructing the evolution operator 

In order to calculate the unitary evolution operator U{t,to) we will perform a generalization 
of the method developed in [2l[ that will clarify the appearance of the auxiliary Ermakov-Pinney 
solution (|2.10p in this context, and will allow us also to warn the reader about other problematic 
choices that have appeared before in the related literature. We first introduce on the (one- 
parameter family of) unitary operators 

D{x) := exp (^-^^jQ^^ , x G M , 

generating a displacement of the momentum operator, D{x)PD~^{x) = P + xQ (the position 
operator being unaffected by them), and define the unitary squeeze operators 



S{y) := exp (^^viGP + PQ)^ , y^M, 



scaling both the position and momentum operator as S{y)QS ^{y) = e^Q and S{y)PS ^{y) = 
e~yp, respectively. Let ^(t) G ^H{t)i t G /, be a solution to the Schrodinger equation, i.e., 
i^{t) = i?(t)^(t), and let x,?/ G C^i^) ■ We now introduce the unitary operators 



T{t)=T{t;x,y) := S{y{t))D{x{t)) , 
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where the functions x and y remain arbitrary at this stage. Let us consider the time evolution for 
the transformed state vector 

«>(t) = $(t;x,y) ■.= T{t-x,y)-^{t), 

given by 

= (T{t)H{t)T-^{t) -iT{t)t{t))^{t) 



= 1 (^e-2j/{t)p2 ^ _ y^t)){QP + PQ) + e2^W(x2(t) + ^(t) + x(t))Q2^ ^>(t) . 

We note at this point that it is possible to get a notable simplification of the previous expression 
just by imposing 

x{t)=y{t) and ^^(t) + + = exp(-4y(t)) . (3.14) 

The most natural way to achieve this is to choose 

y{t) := log p{t) and, hence, x{t) = p{t)/ p{t) , 

with p being any solution to the auxiliary EP equation (12. 9[) introduced in section\l^ In this way, 
the state vector <I>(t; p/p, log /?) =: ^p{t) satisfies the differential equation 

Solving this equation and going back to the original state vector \I'(t), we finally obtain the unitary 
evolution operator for the system. We can then enunciate the following theorem. 

Theorem 3. The time evolution operator U{t,tQ) for the quantum TDHO whose dynamics is 
governed by the Hamiltonian \3. 9^) is given by a composition of unitary operators 

Uit,to) =T-\t)Rpit,to)Tpito) , 

where 

Rpit,to) := exp (-^ / ^(^' + Q')) ' (^.15) 

and Tp{t) = Sp{t)Dp{t), with 

Z)p(t) :=exp(^-^^Q2^ and 5,(t) := exp (^^ logp(t)(QP + PQ)^ (3.16) 

Remark 5. Note that instead of introducing p, we could have used other choices for the x and 
y functions. In these cases, conditions (j3.14p may not hold and the expressions of the evolution 
operator would differ from the one obtained here. For instance, one can select x{t) = u{t)/u{t) 
and y{t) = \ogu{t) as in [21J, with u{t) being any solution to the TDHO equation, but this choice 
is problematic because the set {t G / 1 u{t) = 0} must be non-empty and, hence, the resulting 
formula for the unitary operator is generally not well-defined for all values of the time parameter 
t. This is the reason why the election of the Ermakov-Pinney solution is especially convenient in 
this context -recall that p is a positive definite function. It follows from the above argument that 
the appearance of this solution is nearly unavoidable in this context. 
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Note that the eigenstates of the Rp{t, to) operator (j3.15p are given by those of the Hamiltonian 
operator corresponding to a quantum harmonic oscillator with unit frequency n{t) = 1, 

Ho := ^{P^ + Q^) . (3.17) 

This fact will be shown to be particularly useful to calculate the Feynman propagator. It is also 
important to point out that the procedure employed in this section is implicitly based upon the 
transformation of the so-called Lewis invariant [S] 

Ip{t) ■■= \ + [p{t)P - pmf^ , ipH = 0, (3.18) 

into an explicitly time-independent quantity -although in order to obtain the unitary operator it 
has not been necessary to use it. In this case, we simply have 

T,{t)Ip{t)T^\t) = Ho. (3.19) 

The Lewis invariant is often used to generate exact solutions to the Schrodinger equation, and turns 
out to be especially useful to construct semiclassical states for these systems, as will be discussed 
later. 



C. Propagator formula 

We finally proceed to derive the Feynman propagator for the quantum TDHO corresponding to 
the Hamiltonian (j3.9p . In the previous subsection, we have written down the evolution operator 
for this system explicitly in closed form in terms of the position and momentum operators (see 
theorem\^. It is given by the product of the unitary operators p.l5p and (j3.16p . We calculate now 
the action of these factors on test functions Tp S ^(M) C L^(M, dg) in the standard Schrodinger 
representation. First, it is straightforward to see that 



{Tp{t)i;){q) = Vp(t}exp \^--p{t)p{t)q^ j xP{p{t)q) = J^K+ iq,t; qo)^{qo) dqo , 
{T;\t)i,){q) = -±=e^^(^-P^^q^^^{q/p{t)) = j^K-^ 

where we have introduced the distributions 

K+{q,t;qo) := T^^exp (^-^p(t)p(t)g2^ 5(^0 - pCO?) , (3.20) 
K;{q,t;qo) := -i= exp (^^g^^ %o - ^M*)) • (3.21) 

The propagator for Rp(t,tQ), satisfying 

{Rp{t,to)'ip){q) = / Kp{q,t;qo,to)'ip{qo)dqo , 



can be easily derived from the one corresponding to the TIHO with unit frequency. As is well 
known 33|, |3j] , the Green function for the Hamiltonian (|3.17p is given by the Feynman-Soriau 
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formulae 



K'^{q,v;qo,0) = ^L=sm ^/^(u,0)exp ( — — ({q^ +qQ)cosv -2qqo)] , v^O (modvr), 
V^TTi \2smv\ )) 

I ( % Sin \) \ 
r;; (/O) 0) = cos~"'^"(i;, 0) exp ( — J (5((?o — (?/ cos f ) , u = (modvr), 

y COS X) J 

where the so-called Maslov correction factor [sJ], which allows the calculation of the propagator 
beyond the caustics {u £ M : sin(T;) = 0} = {vr/c : A; G Z} , has been conveniently absorbed into 
the definition of sin^/^(i;,0) and cos^/^(f,0) given in the formulation of theorem\^ In view of Eq. 
(j3.15p . we simply get 

t; qoM) = [q, ^ ; go, . (3.22) 



Therefore, 



where 



{U{t,to)^li){q) = {T;\t)Rp{t,to)Tp{tQ)i?){q) = / K(g, t; go, io)V'(9o) d(?o , 



K{q,t;qo,to) = / K {q,t;q2)K^{q2,t;qi,to)K^ {qi,to;qo) dqidq2 ■ (3.23) 

By combining (fOnjl - ffOHl ) with (l2lTl) and (f2l^ . we find the formula for the propagator (IXTnll 
enunciated in i/ieorem[3 expressed in terms of the c and s solutions to the classical TDHO equations 
(j2.ip . As expected, the propagator -and hence the evolution operator itself- does not depend on the 
particular solution p to the EP equation (j2.9p chosen to factorize C/(t, to)- Taking the appropriate 
limits one obtains, after straightforward calculations, the propagator evaluated at caustics (j3.1ip . 
The resulting expressions are in agreement with those obtained by other authors (see, for example, 
[2^ . [2^ . 34 . HH]), though in our case they have been attained within a different scheme, based 
essentially on the previous obtention of a closed expression for the evolution operator. Finally, a 
direct calculation shows that the propagator K{q,t;qQ,tQ), viewed as a function of {q,t), formally 
satisfies the evolution equation 

idtK = -^d''gK+^q\{t)K. 



D. Transition amplitudes and vacuum instability 

The exact expressions for the Green functions (j3.10p and (j3.1ip can be used to exactly com- 
pute both transition amplitudes and probabilities. Here, we will restrict ourselves to the class of 
normalized states i;^^ defined in L^(M, dg) as 

^1/4 , . 

C(g) := ^==exp -u;gV2U^n(V^(/), u; > , n G No , (3.24) 

with Hn{z) denoting the nth Hermite polynomial in the variable z . Note that for any fixed value co 
the set (0^ : n G No) defines the usual orthonormal basis of L^(M) constituted by the eigenvectors 
of the quantum Hamiltonian ()3.9p corresponding to a TIHO of constant frequency = u. 
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Since the 0!;^ states are complete, the corresponding transition amphtudes and probabihties for 
other states are readily obtainable. By using the generating function for Hermite polynomials, 



exp {2^qx - x^) = ^ , 



n=0 



it is clear that 



(C |C^(i2,ti)C) = - { ''''^''^f^!^^ ^''' s-^^^ (3.25) 



where x^2\f(x\,X2) denotes the complex coefficient appearing in the ^"^Xg^-term of the Taylor 
expansion of the function /. Here, for any matrix A G Mat2x2(C), we define 



l{xx,X2\!y) := exp ( - (x? + X2)) J ^exp (^-^q^Aq + 2x^ diag{^/uJl, ^/LJ^) ^ d^q 
= ^^^^^-^ exp ^x^ (^2diag( A-"Miag( v^) - l) > 



whenever Re (A) > and detA 7^ . In this formula, x denotes the column vector with first and 
second components given by xi and X2, respectively; we define q similarly. In our case. 



/ .c(t2,ii 

LUi — I 



A{ti,t2;uJi,U2) := 



V Sit2,ti) 



S{t2,ti) s{t2,ti) 

i .C(ti,t2 
i02 - i 



with 



det A(ii,t2;<^i,W2) = (^iuJ2 - ~ ^ ^ 



s{t2M) / 

UJlc{ti,t2) +UJ2C{t2,ti) 
S{t2,tl) 



Here, Re{A{to,t;oji,oj2)) > and det A(to,t;wi,W2) 7^ for all {to,t) e Ixl anda;i,ti;2 G (0, +00) . 
The Taylor expansion of I{xi,X2; A) can be efficiently computed by applying the following lemma, 
that trivially follows from the multinomial formula. 



Lemma 1. Let 



B = B* 



m fi2 

'12 ^'22 



G Mat2x2(C). 



Then, using the notation introduced above, we have 

exp =6S7~"^)/^ (2612)"^ V „ , , , 

2^ / , m!(m+(ni-n2)/2)!(n2-2TO)! 



(6ll622)"^(46?2)~'" 



meA(ni,n2) 



whenever Hi and 77,2 have the same parity, and vanishes otherwise. Here, A(ni,n2) := (m G No : 
max{0, (n2 — ni)/2} < m < [^2/2] ) , where \_x\ denotes the largest integer less than or equal to 
X G M . In particular, taking ni = , we get 



[x^^x'i^] exp { x'Bx ) = 



,n2/2 
"22 

W2)! 



for ^2 = (mod 2) , 



(3.26) 
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and vanishes if n2 is an odd number. 

Remarks. Note that the TDHO quantum dynamics is invariant under parity inversion P and 
the states c/)^ satisfy Pc/)^ = Hence, | [/(tg, = if ni and n2 have different 

parity. 

As a concrete example, in the case of a TIHO with constant frequency uj = uji =1^2, 
we identify 



B = 2diag{^/uj, ^/uJ) A ^{ti,t2;uj,uj) diag{^/uj, ^/uj) - I = exp ( - zu;(t2 - ^i)) 
and, hence, 

I(xi,X2;A(ti,t2;uj,uj)) = y — rexp ( — iu;n(t2 — ti))x?X2 . 



1 

1 



n! 

n=0 



This is in perfect agreement with 

;7(t2, ti)Ci> = exp ( - iuj{ni + l/2)(t2 - h))6{ni,n2) , 



"2 



where (5(ni, 71.2) denotes the Kronecker delta. For arbitrary time-dependent frequencies the formula 
(j3.25p . when restricted to the same initial and final frequencies lvi = 002, coincides with the one 
given in [s^ written in terms of associated Legendre functions. 

We conclude this section with the analysis of the instability of the vacuum state (pQ due to the 
nonautonomous nature of the Hamiltonian (13. 9[) . This can be easily derived from the formulae 
(I3:25]l and I^M)- 

Theorem 4. The quantum time evolution of the vacuum state (pQ is generally given by a su- 
perposition of states C/(t,to)0o ~ SneNo^'^an I ^(*' *o)</'o)'?^2n ' where the probability amplitudes 
{<p2n \U{t,to)(p'^) are given by 

{^2n I U{t,to)<P^) = y^^{2Lo{A-\to,t;u;,uj))22 - | U{t,to)cP^) , n G N, (3.27) 

in terms of the the expectation value 



{r, I U{t, toW^) = J r exp(-i7r/4) 5-^2 (t, to) 

y aet A[tQ,t;uj,io) 



with 

.2 



(A-^(t0 t-LO UJ))22 = - isit2,ti)c{t2,ti) 

1+Uj'^s'^{t2,ti) - C{t2,ti)c{ti,t2) - iu;s{t2,ti){c{t2,ti) + C{ti,t2)) 

Remarks. Consider the usual annihilation and creation operators 

:= ^ (V^g + iP/ V^) and a* := ^ (V^Q - iP/V^) , (3.28) 
V 2 V 2 

with [auj,alj] = 1 and [a<^,a^] = = [a2j,alj], such that a* = ^yn + Icp^+i and a^c/)!;^ = ^/n(j)'^_-^^, 
Vn G N, with a^^cpQ = 0. The evolution of these operators in the Heisenberg picture can be 
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obtained directly from Eq. (j3.6p . 

U-\t,to)a^U{t,to) = A^{t,to)a^ + B^{tM)al , (3.29) 
U-^{t,to)alU{t,tQ) = B^(t,to)a^ + A^{t,to)al , 

where A^{t,tQ) and B^{t,tQ) are the Bogohubov coefficients 

A^{t,to) := ]^(^c{t,to)+s{t,to)+i{LO-^c{t,to)-uJs{t,to))y (3.30) 

B^{t,tQ) := ^(c(t,to)-s(t,to)+«(w~^c(t,to)+'^s(t,to))) , (3.31) 

satisfying A^{tM) = A^{to,t), B^{t,to) = -B^{to,t), and \ A^{t,to)\'^ - \B^{t,to)\'^ = 1, V(t,to) e 
I X I . Note, in particular, that Ai^{t,to) never vanishes. For example, for the TIHO of constant 
frequency a; > we have B^{t,tQ) = and A^{t,tQ) = exp(— i(t — tQ)uj). A straightforward 
calculation yields (see also [sjI and [s^ ) 



Uit,to)^^ = (0^ I UitM^^) exp [-\j^^f*') <Po , (3-32) 

This formula is in perfect agreement with the transitions (j3.27p . Indeed, it is straightforward to 
check that 

2Lo{A-\to,t;uj,io))22-l = -BUto,t)/AUto,t). 

Since det A{tQ,t;io,uj) = —2i(jjs~^{t,tQ)Ai^{to,t), the expectation value (^g I f^(^' *o)'Ao ) can be 
rewritten as 

(c/.^ I U{t, to)0^) = ,, ^ \ . , exp(icT(t, to)) , (3.33) 
y'\Auj[to,tj\ 

where the phase a{t,to) G C^(/ x /) comes from a careful calculation of the principal argument. 
For a TIHO with constant frequency w > 0, we have a{t, to) = (to — t)uj/2 for all t, to G IK • Given 
an arbitrary squared frequency K(t), the phase cj(t,to) evaluated at times t close to to is simply 
given by 

. s 1 , /u;s(t,to) - a;-ic(t,to)\ 
ait, to) = -- arctan ( ^(t, to) + i(t, to) J ' ^'''"^ 

The C7 phase can be conveniently canceled through a suitable redefinition of the Hamiltonian 
(13. 9p just in the case when &(t, to) is independent of to . In that situation, by identifying ■i?(t) = 
(j(t, to) in Eq. (|3.12p . we have that the redefined evolution operator satisfies {(pQ \ ?7i(t, to)0o) = 
l/y/\A^{to,t)\ . In the TIHO case, we get 'd{t) = —uj/2 (this amounts to considering normal order). 
In general, it is not possible to proceed in this way in all situations when dealing with arbitrary 
time-dependent frequencies. In any case, the a phase is irrelevant for the calculation of transition 
probabilities. In particular, given ^i, ^'2 G with ^'i = Fi{auj, a* ) (p^ , where Fi is some analytic 
function, we have 

Prob(^'2,t2|^i,ti) = K^'2|f/(i2,ti)^i>P 

^ I (^2 I Fi{a^Hih,t2),aZHiti,t2)) exp ( - ^..(ti, t2)/(2A^(ti, t2))o*^)(/>g)P 

|^c^(tl,t2)| 
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where the time dependence only appears through the Bogohubov coefficients ()3.30p and (j3.3ip . 
Finally, it is important to point out that the transformations ()3.29p and the evolution of the 
vacuum state ()3.32p fully characterize the quantum time evolution of the TDHO. By using these 
relations, we can easily compute the action of U{t,to) on any basic vector 0^!^ = (l/\/nI) a^4>Q . 



IV. SEMICLASSICAL STATES 



In this section, we will look for states that behave semiclassically under the dynamics defined 
by the quantum Hamiltonian (13. 9p . We will base our study on the concrete factorization of the 
evolution operator defined in theorem\^ To achieve this goal, note that the eigenvalue problem for 
the Lewis invariant ()3.18p can be exactly solved. Indeed, let us fix ^ ^-nd let : n £ No) 
be the eigenstates ()3.24p of the auxiliary Hamiltonian Hq (j3.17p corresponding to unit frequency 
(jj = I. According to the relation (j3.19p . the initial states ipn{to) := T~^{to)cl)n, 

1 fi (m 



labeled both by p and the integers n G No, are eigenstates of Ip{to) with eigenvalues equal to 
n + 1/2. Consider now the initial pure state 

oo „ 

<^(;\to) := T-i(to)$(^) = e-l^lV2 ^ A_v,P(to) , zee, (4.1) 

n=0 

with := ■^^^(/'n being the well-known coherent states for the Hamiltonian Hq . 

Let us take the annihilation and creation operators a and a* for unit frequency uj = 1 defined in Eq. 
(j3.28p . The superposition (|4.ip is a normalized eigenvector of the (time-dependent) annihilation 
operator 

apito) := T~\to) a Tp{to) = (Q/p(to) + i{p{to)P - pi^M) , (4-2) 

in the sense that ap(to)^p^\*o) = z^p^\tQ) . This operator, together with the associated creation 
operator 

a;{to) := T-\to) a* r,(to) = ^ {Q/p{to) - i{p{to)P - p{h)Q)) , 

satisfies the Heisenberg algebra, [ap(to), a*(to)] = 1 and [ap{to), ap{to)] = = [a*{tQ),a*{tQ)], for 
each initial value of the time parameter to • In particular, the Lewis invariant (j3.18p may be 
expressed in terms of these operators as /p(to) = '^p(^o)op(to) + (1/2)1- Through unitary time 
evolution, we get 

$(^)(t,to) := U{t,to)<^'^;\to) = exp (-^ jj" ci>(^''(*'*«»(t) , (4.3) 



where we have denoted 



Zp{t,to) ■■= exp (^-i p^^) ^' ^^'^ 
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We want to remark that the time-dependent phase appearing in Eq. (j4.3p is necessary for these 
states to verify the Schrodinger equation. In our case, they coincide with those defined in Eq. 
(4.6) of reference [sof]. We conclude that the family of states ()4.ip is closed under the dynamics. 
Moreover, the following theorem can be used to justify that these states can be considered as 
semiclassical under certain assumptions. 

Theorem 5. Let z = x + iy £ C and Iq £ I. The position and momentum expectation values in 
the state ^l,''\t,to) = U{t,to)^l,''\to) satisfy 

qH{t,to) = ($(")(t,to) |Q^i,'^(t,to)> = V2/9(t)Re(zp(t,to)) , 
PHit,to) = {<S>'-;\t,to)\ P<^>'-;Ht,to)) = V2Re{{pit) - i/p{t))zpit,to)) , 

where (qhtPh) is the classical solution iS. 2\) determined by the Cauchy data {q,p) = 
[V2p{tQ)x, \'^{p{tQ)x + y / p{to))) at time to . Moreover, the mean square deviations of the position 

and momentum operators with respect to the evolved state $p (t, io); 

V'(*,*o)^ = ^ Pit) , \i^\tMP = 71 - 'P''i')\ ' (4.4) 
are independent of both to and the Cauchy data defined by z . 

Remark 6. Given any observable O, its uncertainty in the state ^ G QIq is defined as AipO := 

1 /2 

((^' I O^^') — I O^')^) . Note that, in general, the elements of the family of states under 
consideration are neither standard coherent states nor squeezed states. For instance, for the free 
particle (12. 4p one can choose p{t) = y^l + (t — to)^ and, hence, A^(z)^^ ^^^Q ~ t/\/2 for large values 
of t ; similar results occur for other elections of p . Nevertheless, it is obvious that we will obtain 
good semiclassical states for a system whenever the solution p to the auxiliary EP equation (12. 9p 
has a suitable behavior, for instance, if p is periodic in time or is simply a bounded function. We 
will analyze some clarifying examples in this respect. 

Example 1 (Vertically driven pendulum). Consider the vertically driven pendulum [4ol |. i.e., the 
motion of a physical pendulum whose supporting point oscillates in the vertical direction. In the 
small angles regime, it is described by the Mathieu equation in its canonical form [41,] 

'u(i) + K;(t; a, 6)'u(i) = , K(t; a, 6) := a - 26cos(2t) , a, 6gM. 

The general solution to this equation is a real linear combination of the so-called Mathieu cosine 
and sine functions [i^, S^, denoted respectively as Ce{t;a,b) and Se{t;a,b). Given a nonzero 
b value, it is a well-known fact that the Mathieu cosine and sine functions are periodic in the 
time parameter t only for certain (countable number of) values of the a parameter, called char- 
acteristic values. The procedure to calculate these characteristic values for even or odd Mathieu 
functions with characteristic exponent^ r € Z and parameter b can be efficiently implemented in 
a computer. In this case, solutions to the EP equation (j2.9p inherit the periodic behavior from 
the Mathieu solutions, in such a way that one obtains well-behaved semiclassical states for which 
the average position and momentum follow the classical trajectories, whereas the corresponding 
uncertainties vary periodically in time. Note that, for small values of the b parameter, we have 
Ce{t;a,b) ~ cos(-y/at) and Se{t;a,b) ~ sm{^/at), and the system closely approximates the TIHO 



^ All Mathieu functions have the form exp{irt)F{t), where r is the characteristic exponent and the function F{t) 
has period 2n. 
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FIG. 1: Variances of the position and momentum operators for the 3-torus Gowdy-type oscillator. Here, 
p{t) = ^ irti^jQiujt) + Y^{ujt)^ /2. The ^\^\t, to) are states of minimum uncertainty for times t far from the 
singularity at t = . 



with squared frequency given by the a parameter. 

Example 2 (T^ Gowdy-like oscillator). Consider the TDHO equation 

u(t) + K(t;uj)u(t) =0, K(t;uj) := uj'^ + ^ , u; € M , t € (0, +oo) . 

This equation is satisfied for each mode of the scalar fields encoding the information about the grav- 
itational local degrees of freedom of the so-called Gowdy models, which are symmetry reductions 
of general relativity with cosmological interpretation that admit an exact -i.e., nonperturbative- 



quantization (see next section). In terms of the zero Bessel functions of first and second kind 43l |. 
denoted Jq and Yq respectively, the c and s solutions introduced in section O are given by 

c{t,to) = ^U ^Yo{u;to)-2ujVt^Yi{u;to)yUujt) - ^ ( JqMo) - Sc^^ti^JiMo)) i^oM), 



s{t, to) = -^VhiYo{Loto)Jo{ujt) + ^Vt^Jo{uJto)Yo{ujt) . (4.5) 

Note that the squared frequency is a sum of a positive constant plus a decreasing function of 
time, so that the system approaches a time-independent oscillator as t tends to infinity. In Fig. 
[H we show states $p (t,to) that behave as coherent states for large values of the time parameter. 
The classical equation of motion has a singularity at t = which translates into the vanishing of 
the uncertainty of the position operator -and, hence, into the divergence of the variance for the 
conjugate momentum- at that instant of time. 

There are other interesting effects due to the classical singularity. Let us consider again the 
study of transition amplitudes developed in subsection \IIID\ and take uji = uj2 = OJ. We proceed 
to analyze the behavior of the (unique) state \E'(f2,ti) that evolves to the vacuum state (p^ at time 
t2 when used as Cauchy data in ti <t2 , i.e., 

U{t2M)^{t2M) =^0^ *(t2,tl) = C/(tl,t2)</'0 ■ 

The transition amplitudes (02n I ^(^2,^1)) = (02n I ^(*i)*2)0o)) ^ I^O; can be computed by using 
Eq. (j3.27p . We recognize two regions of interest in the time domain, 

To+:={{ti,t2)\0<ti<.uJ-^ <.t2} and r++ := {(ti, ^2) | ^^"^ < ii < ^2} • 
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0' ' ' t I ' ' t 

^ n 1 ^ 

FIG. 2: Variances of the position and momentum operators for the 3-handle and 3-sphere Gowdy-type 
oscillators. Here, we take the solution p{t) = \J sin t ( S^f^^, _ i) /2 (cos t) + =S^t^'_i)/2(cost))^^^ to the auxiliar 
Ermakov-Pinney equation. In particular, these graphics correspond to cj' = 5 . 



In r++, the asymptotic behavior of the Bessel functions for large values of the time parameter 
[4^ leads the system to behave as a TIHO of constant frequency a;, with ^'(t2,ti) ~ . On the 
other hand, in the region Tq^ , the proximity of ti to the classical singularity manifests itself in 
the fact that the wave function takes the form ^{t2,ti) ~ . Note that this behavior is in conflict 
with the unitary evolution of the system, which implies ||^'(t2, ti)|| = 1 • 

Example 3 (S^ x and S'^ Gowdy-like oscillators). Gowdy models admit spatial topolo- 
gies different from the 3-torus one, concretely the 3-handle x and the 3-sphere S^. As 
expected for closed universes, these systems present both initial and final singularities. For this 
reason, they become useful test beds to discuss the exact quantization of cyclic universes. Here, 
the modes satisfy equations of motion of the form 

u{t) + K{t;uj)u{t) = 0, Hi{t;uj) ■.= uj'^ + ^{l+csc'^t) , cjGM, t £ {0,tt) . 

In this case, in terms of first and second class Legendre functions ^43i] denoted respectively as 
and , X G M , we have 



c{t,to) = - \J sin 1 1 sin to^(a;'-l) /2 (cos t) ( (1 + w')=^(l+a;')/2(cOSto) " t^'^{a;'-l)/2 (cOS to)) 

- ^ sin tj sin to d^{uj' _ i)/2 (cos t ) ( ( 1 + w')^(i+a;')/2(cos to) - w'^(^/_i)/2(costo)) , 

s(t,to) = Vsintsinto(^(a;'-l)/2(cOSto)=^(a;'-l)/2(cOSt) - =^(a;'-l)/2(cOSto)^(a;'-l)/2(cOSt)) , 

(4.6) 

where u' := \/l + 4a;^ . In Fig. [21 we show the behavior of states <I>p^^(t,to) for which the uncer- 
tainties of the position and momentum operators have an oscillatory behavior far enough from the 
singularities occurring at t = and t = vr. Although p does not vary periodically, the function 
remains bounded and, thus, the <I>p {t,to) states can be used to perform a semiclassical study of 
these models. Finally, one may proceed as in the 3-torus case in order to analyze the way the 
classical singularities affect the quantum behavior of the systems, obtaining similar results. 
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V. EXTENSION TO FIELD THEORIES 
A. General framework 

In this section, we will extend the previous study to linear dynamical systems with infinite de- 
grees of freedom, focusing our attention on the study of scalar fields evolving in fixed backgrounds 



such as the Minkowskian space or the so-called (linearly polarized) Gowdy metrics [2^, that cor- 
respond to symmetry reductions of general relativity describing cosmological models with initial 
and final singularities.^ Our results, however, will have a wide range of applicability, being pos- 
sible to easily extend them to other field theories with similar structure. The construction of the 
appropriate C*-algebra of quantum observables can be obtained by making a simple comparison 
with the one-dimensional case discussed in section ing Consider the canonical phase space (r,a;) 
consisting of the infinite-dimensional M- vector space T of smooth Cauchy data endowed with the 
natural (weakly) symplectic structure u. Taking the linearity of F into account, each element A G F 
is identified with the functional : F — > M such that, for any other A' € F, Fx{X') := cj(A, A') . 
The abstract quantum algebra of observables is then given by the Weyl C*-algebra on F, #'(F), 
generated by the elements W{X) = exp(iFx), A G F, that formally satisfy the relations ()3.3p . The 
GNS construction [l3] establishes that, given any state loq : W{T) — > C on the algebra -that is, 
a normalized positive linear functional-, there exist a Hilbert space {J%, {'\')jifb)^ ^ representation 
ttq : W{r) — > ^(J^) from the Weyl algebra to the collection of bounded linear operators on J%, 
and a cyclic vector £ =^ such that 

uJo{A) = (v&o I ^o(^)^o),if& , e W{T) . 

Moreover, the triplet {J%,ttq,"^q) satisfying these properties is unique up to unitary equivalence. 
Von Neumann's uniqueness theorem, however, cannot be generalized to field theories, and states 
on the Weyl algebra generally yield nonequivalent cyclic representations. In order to pick out a 
preferred representation of the canonical commutation relations one can impose additional crite- 
rions, such as the well-known Poincare invariance of the Fock representation for scalar fields in 
Minkowskian space. Regarding the Gowdy models, one imposes the invariance under an extra 



U{1) symmetry generated by a residual global constraint for the 3-torus case 4j], or the invariance 



under the symmet ry g roup 50(3) of the Klein-Gordon equations of motion for the 3-handle and 



the 3-sphere cases |27I]. 

Every linear symplectic transformation T G SP(F), for which a;(TAi,TA2) = a;(Ai,A2), 
Ai, A2 € F , defines a unique *-automorphism ar G Aut(^(F)) such that (ar o W){X) := T^(TA) . 
This is the case, in particular, of the symplectic transformations that characterize the classical 
dynamics of the system. Given a concrete Hilbert space representation (.y^,7ro,^o) of the Weyl 
C*-algebra, the symplectic transformation T G SP(F) is said to be unitarily implementable on the 
cyclic representation space J^o if ttq and nooar are unitarily equivalent, i.e., there exists a unitary 
operator Ur '■ J% such that 

U^^7t{W{X))Ut =TT{{aT oW){X)) , yW{X) eW{r). (5.1) 

A common feature of the quantization of infinite-dimensional linear symplectic dynamical systems 
is, precisely, the impossibility of defining the unitary quantum counterpart of all linear symplectic 



^ The isometry group of the Gowdy models is 17(1) x (7(1) and the spatial slices are restricted to have the topology 
of a 3-torus T^, a 3-handle x S^, a the 3-sphere S^, or the lens spaces L{p,q) (that can be studied by imposing 
discrete symmetries on the 3-sphere case). The exact quantization of the linearly polarized Gowdy models in the 
vacuum or coupled to massless scalar fields has been profusely analyzed (see [2^ and [S^l and references therein). 
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transformations on the phase space [45[]. This is the case of the time evolution of the Gowdy 



models when these systems are written in terms of the dynamical variables that naturally appear 
after performing their Hamiltonian formalisms. Note, however, that the lack of a unitary operator 
implementing the quantum time evolution conflicts with the axiomatic structure of quantum theory 
itself. In case of not rejecting the models for this reason, one must carefully analyze the viability 
of a suitable probabilistic interpretation for them, as discussed in [4^. Nevertheless, it is possible 
to overcome this problem by performing some suitable time-dependent redefinitions of the basic 
fields [25|, l27|. In what follows, we will always refer to these rescaled fields. 



B. Unitary quantum time evolution 

The canonical phase space description of the classical systems under consideration consists now 
of an infinite-dimensional nonautonomous Hamiltonian system (/xF, dt, u, H(t)). Here, T := '^y.^ 
is the space of Cauchy data, where denotes the Frechet space of rapidly decreasing real sequences 



X := {xi : £ G X), with i running over some countable set X 47|, l48l . l49l |. This space is endowed 
with the natural symplectic structure 

^((qi,Pi), (q2,P2)) := XI ip^^it - V2l(l\l) , V(qi,pi), (q2,P2) G T. 



The time-dependent Hamiltonian H{t) : F ^ M is a quadratic form on F that can be diagonalized 
as a sum of TDHO Hamiltonians of the type ()3.ip . Explicitly, 



H{t,q,p):=^^{pj + K,{t)qj), i G / = (t_ , t+) C M , (5.2) 



2 



where the time-dependent squared frequencies K£{t) G C^{I), £ G X, must satisfy K{t) = {^^{t) : 
£ G X) G for all t €z I, with ^ denoting the vector space of slowly increasing real sequences, 
i.e., the topological dual of The time evolution is implemented by symplectic transformations 
of the type ([321) 



where, for each ^ G X, q and si are the solutions to the TDHO equation of squared frequency ne,{t) 
introduced in section\ll\ 

For the Minkowskian quantum field theory generalized to a spacetime R x T'^ with closed spa- 
tial sections we have X = Z \ {0} and ^^(t) = for all t G M. For the 3-torus Gowdy models 



25| we take X = Z \ {0} and Ki,{t) = + with t G (0,+oo). For the remaining 

topologies admitted by the Gowdy cosmologies, the 3-handle and the 3-sphere [23], X = N and 
Ki{t) = + 1) + (1 + csc^ t) /4 , with t G (0, vr) . For simplicity, we will not consider in these cases 
the quantization of the zero mode ^ = 0. It can be represented in terms of standard position and 
momentum operators with dense domains in L^(M) . The q and functions for the Minkowskian 
and Gowdy 3-torus cases are respectively given by (j2.3p and (j4.5j) substituting = ; for the 
3-handle and 3-sphere topologies, they are given by (j4.6p identifying uj = \l i{t + 1) . The classical 
singularities of the Gowdy models will persist in their quantum formulations. 

In order to exactly quantize the infinite-dimensional systems under consideration, we introduce 
Schrodinger representations S^, 51) 52, [i^, where state vectors act as functional ^' : ^ ^ C be- 



longing to certain Hilbert spaces M'ol = L {"^^ o"(Cyl(^)), d/Uct) ) and define suitable decompositions 
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of the position and momentum operators in terms of modes. There are subtleties associated with the 
infinite-dimensionahty of the classical configuration space that affect the definition of the Hilbert 
space. On one hand, it is not possible to define nontrivial Lebesgue-type translation invariant mea- 



sures /ict, but rather probability ones |49l. l54i| . On the other hand, the set over which the measure 
space is built -the so-called quantum configuration space- must be given by some suitable distribu- 
tional extension of . In this case, it suffices to consider the dual The reason to proceed in this 
way is that the measure is not supported on i.e., the classical configuration space has zero mea- 
sure with respect to iia.. Given a nonzero complex sequence cx = (a^ : £ € X), the Gaussian mea- 
sure fj,a is defined on the cylinder sets cr-algebra c7(Cyl(^)) on ^. This is the smallest cr-algebra with 
respect to which the functionals q S ^ i— > (q, x) := X^^gx H^t- measurable for each x G Ex- 
plicitly, for each finite n-tuple {£i, . . . , in) G X", such that ii < {i = 1, . . . , n—1), let us consider 
the projections pi^-e^ ■ — > M", q i-^ pi^...i^{q) = (qi-^, . . . ,qi^) ■ Then, fia is defined by its action 
on cylinder sets belonging to u-algebras of the form = pj^...^ (-B(M")) C (j(Cyl(^)), 

where B(W^) is the Borel c-algebra of subsets of M", i.e.. 

The sequence a defines a positive continuous nondegenerate bilinear form Ca : ^ x — > M through 
the formula Cci(x, y) := Xl^ex I'^^P^^^^ ' x, y S ^. Ca is called the covariance operator in this 



context [481]. 



We consider now certain class of measurable functions that we will use in the following. Let 
C ^ be a finite-dimensional subspace of the classical configuration space A cylinder function 
^, based on is a function of the form 

^'(q) =F((q,Xi),...,(q,x„)) 

for a finite set xi, . . . ,x„ G J^, where F : M" — > C is a smooth function. They are called cylinder 
because they depend on q G through the pairings (q,Xj) defined by a finite number of "probes" 
Xj G ^. The Gaussian measure fia can be used to endow the linear space spanned by this type of 
functions with an inner product. Explicitly, in the particular case ^ = spanje^^, . . . ,e^^}, where 
ej G ^ denotes a sequence whose only nonzero component is the ith one, every cylinder function 
^' on ^ can be written in the form ^'(q) = F(qi-^ , ■ ■ ■ , qe„) and the scalar product is given by 



„ n 

(^'l|^'2)= / Fi{qe„...,qijF2{qi„...,qejYl 



1 

exp 



2tt \ ai. 




Square integrable cylinder functions on ^ span a C- vector space that we denote as Cylj;-. The 
class of all cylinder functions is denoted by Cyl := U^Cyl^. The inner products defined on each 
Cyl^ can be extended to Cyl in the natural way and the Cauchy completion of Cyl with respect 
to this inner product is the Hilbert space Jifa (see [si*] for more details). 

The configuration observables will act as multiplication operators, whereas the canonically con- 
jugate momenta will differ from the usual ones in terms of derivatives by the appearance of multi- 
plicative terms which are necessary to ensure their self-adjointness; specifically, for ^ G Cyl, 

(Q,M/)(q)=(?,^(q), (P,xl/)(q) = -i^(q) + |^g,.l/(q). (5.3) 

oqe ai 
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The complex sequences a and (5 must satisfy 

a^Pi-P^ai = i, WieX, (5.4) 

by virtue of the CCR, [Qe,Pei] = i6{£,£')l and [Qi,Qe'] = = [P£,P£/]. These conditions imply 
lo^ll/^^l > 1/2 for all ^ G X. The functional form of a and (3 for the Minkowskian case and Gowdy 
models will be discussed later in the context of the unitary implementability of their dynamics. 

According to condition (15. ip . if the quantum dynamics is unitarily implementable there exists a 
(biparametric) family of unitary operators C/(t, to) : -^a, depending on (t, to) & I ^ I, such 

that 

U-^{tM)QiU{t,to) = Ci{t,tQ)Qe + stit^Pt , (5-5) 

The above relations characterize to) univocally up to phase. They can be rewritten in terms of 
annihilation and creation operators and a| , with [a^, a|,] = 5{l, and [ai, a^/] = = [a|, a|,] , 
such that 

Qi = aiae + a^a^ , Pi = Peae + Pea*e <^ ag = -iPeQe + iaePe , = il3iQe - iaePe ■ (5.7) 

Relations ()5.5p and (|5.6p are then equivalent to 

U^\t,to)aeU{t,to) = Ai{t,to)ae + Be{t,to)a} , (5.8) 
U-\t,to)a*eUit,to) = Beit,to)ai + Ai{t,to)a}, 

where the Bogoliubov coefficients Ai{t,to) and B£{t,tQ) are given by 

Ai{t,to) := i(^Si{t,to)aiPi - Ci{t,to)(3iai + ce{t,to)\ae\'^ - Si{t,to)\Pi\'^^ , (5.9) 
Bi{t,to) := i(^{se{t,to) - Ci{t,to))a£P£ + C£{t,to)aj - Si{t,to)P'^^ , (5.10) 

satisfying \Ai{t,to)\'^ -\Bi{t,to)\'^ = 1, for all £ G X . In particular, Ai{t,to) / for all (t,to) & Ixl ■ 
According to the theory of unitary implementation of symplectic transformations [IH] , the unitary 
time evolution operator [/(t, to) exists if and only if 

B(t, to) = (Beit, to) : i\C) ^ Yl *o)l' < +oo , V (t, to) G / x I . (5.11) 

For quantum free fields in Minkowskian spacetime the usual choice 

1 f\e\ 

ai = —= and /3^ = -iW-LL (5.12) 

implies B(t, to) = and, as is well known, the time evolution is unitarily implementable for inertial 
observers. The uniqueness of this representation can be proved by using the same techniques 
employed in the standard M x Minkowskian spacetime under the condition of Poincare invariance. 
Concerning the Gowdy models, it is straightforward to show that the unitarity of th e q uantum 
time evolution is guaranteed for a and /3 sequences with the asymptotic expansions (see [25|,l26|,l27| 
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for more details) 



a, = —=exp(i7,) + 0(1^1-3/2) , = /U exp(i7,) + 0(K|-V2) , 
2,\l\ V I 

where 7 is an arbitrary real-valued sequence. Moreover, all {U{\)- or 50(3)-invariant) representa- 
tions for which the dynamics is well-defined and unitary are unitarily equivalent. In what follows, 
we will assume the use of the particular choice of a and (3 coincident with Eq. (j5.12p . As a 
counterexample to the previous cases with unitary evolution, consider a system of infinitely many 
harmonic oscillators with imaginary frequency of the type (|2.5|) . In this case, it is possible to 
show the nonexistence of sequences a and /3 satisfying (j5.4p such that the dynamics is unitarily 
implemented. The 'wrong sign' of the squared frequency is also responsable for the failure of the 
unitarity of the time evolution in more com plic ated systems, such as minimally coupled massless 
scalar fields evolving in de Sitter spacetime [551]. 

Note that the annihilation operators are given by the derivatives 

d 

dqe 

and, hence, the vacuum state ^'o G J^a satisfying a^^'o = for all ^ E X is given by the unit constant 
functional ^'o(q) = 1 up to multiplicative phase. The states with finite number of particles, which 
are obtained as the image of any polynomial in the creation operators acting on the vacuum state, 
define a common, invariant, dense domain of analytic vectors for the configuration and momentum 
operators (15. 3p . so that their essential self-adjointness is guaranteed and, hence, the existence of 
unique self-adjoint extensions (see Nelson's analytic vector theorem in [5fl]). 

From the Bogoliubov transformations (j5.8p . and proceeding as in the one-dimensional case 
(j3.32p . one easily computes the evolution of the vacuum state, 

U{t,to)^o = {'fo\U{t,to)^o)exp (-i^l^jMaf ) , (5.13) 



with 



/ex V\Mto,t)\ 

Due to the unitary implementability of the dynamics, the square summability of the sequence 
{B£{tQ,t)/Ai{to,t) : £ £ X) and the convergence of X^^gx ^^S l^^(*0) *)| are guaranteed and, hence, 
the action of U{t,to) is well defined over states with finite number of particles. The phase of the 
expectation value (j5.14p . though being irrelevant to answer most of the physical questions, can be 
explicitly calculated once a quantum Hamiltonian has been fixed. The Hamiltonian verifies Eq. 
(j3.8p and coincides with the operator directly promoted from the classical expression (j5.2p modulo 
an arbitrary t-dependent real term proportional to the identity which encodes the choice of U{t, to), 
i.e., 



2 



where the sequence '&i{t) G C°(I), £ G X , is usually employed to avoid the appearance of infinite 
phases. Analogously to the one-dimensional case, when the dynamics is unitarily implementable. 
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we define the time evolution propagator through the relation 

(?7(t,to)*)(q) = /_K«/3(q,t;qo,io)^(qo)dAia(qo), 
where a straightforward calculation formally provides 

-?^«^(q,i;qo,io) = \\'^\oti\ey:.v[i:[—qlt-^qi\]Ki,[q^ f dT'diir) 

\^ '^l J J \ Jto 

(5.15) 

with Ki denoting the propagator (j3.10p associated with the one-dimensional oscillator of squared 
frequency ne{t). The reader may wish to compare this expression with Eq. ()3.13p corresponding to 
a single oscillator. This formula coincides with the one obtained in ^] when restricted to generic 
finite-dimensional and time-dependent linear Hamiltonian systems. For the sequences ()5.12p . Eq. 
()5.15p provides the expectation value 



with 



/ex V\Mto,t)\ V ^ Jto 



, ^ 1 , f\£\se{t,to)-\e\-^ce{t,to) 
ciit, to) = — arctan 



Mr) 



2 V ce{t,to) + se{t,to) 

for times t close to to- For the Minkowskian free fields, as in the one-dimensional case, the phases 
ae{t,tQ) can be exactly canceled for all £ G Z \ {0} just by defining a normal ordered quantum 
Hamiltonian, which amounts to choosing ^^^(t) = —\£\/2. In general, however, it is not possible 
to eliminate them. For the Gowdy cosmologies, attending to the fact that the q and si functions 
tend to those corresponding to the Minkowskian case for large values of \£\, it is easy to check that 
normal ordering allows only the cancelation of the phases at high frequencies, with 'd£{t) ~ —\i\/2 
as \£\ — > +00 . 

Since we are dealing with systems of infinite number of uncoupled oscillators, one would expect 
that the analysis developed in section for a single oscillator would allow us to factorize the 
evolution operator in the form 

U{t,to) = T-\t)Rp{t,to)Tp{to) , (5.16) 

where, given a sequence p{t) = {pi{t) : £ G X) of solutions to the auxiliary Ermakov-Pinney 
equations pi + Ki{t)pe = l/pf , the Tp{t) and Rp{t,to) operators are univocally characterized up to 
phases by their action on annihilation and creation operators, 

Tp^{t) aiTp{t) = i (^PittiPiit) - - \ai\'^pi{t)^ ai + i {^(jii{pi{t) - -J^) " > 
Rp^{t,to)aiRp{t,to) = ^cos pf(r)) ~ ^^^'"^'^ + '^^'^^ P^m)) 
- i{aj + /3|) sin ( [ -^rr] 4 > 



.2 

to ■ 



and similarly for a| . Again, the resulting unitary evolution operator should be independent of 
the particular choice of p since, as we have shown above, the propagator does not depend on it. 
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However, even in the case of U{t, to) being well-defined as unitary operator, the factorization (j5.16p 
may be ill-defined. This is, in fact, the case for free fields evolving in Minkowskian and Gowdy- 
type spacetimes, as we will prove below. Obviously, this does not prevent us from defining other 
well-defined factorizations for C/(t,to) different from (|5.16p . A particularly convenient choice is the 



one made in j26| for the 3-torus Gowdy model. Calculations developed there can be translated 
essentially unchanged into the remaining Gowdy spatial topologies and the Minkowskian case 
(see the appendix [All . Here, however, we are interested in the original factorization due to its 



implications for the search of semiclassical states. According to 45|], the necessary and sufficient 
condition for Tp(t) to be unitary for each value of t is given by 

J2\(XiMP£it)-'^/peit)) -ajpeit)]^ <+(x>, VtG/. (5.17) 

Similarly, it is straightforward to show that Rp{t, to) is unitarily implementable if and only if 

y (a^ + /3|) sin / ^— <+oo, y {t,to) £ I x I . (5.18) 

Equations (|5.12p for quantum free fields in Minkowskian and Gowdy-type spacetimes lead us to 
conclude that conditions (|5.17p and (j5.18p are not satisfied and, hence, neither Tp{t) nor i?p(t, tg) 
are unitary in those systems. In the case of Rp{t,tQ), this conclusion follows readily, irrespective 
of p{t). For Tp{t), a necessary condition for (j5.17p to be satisfied is given by 

Y^\pi{t) - 1/ pi{t)\'^ <+oo ^ \im pe{t) = l, VtGl, 

where we have taken into account the fact that the real sequence {pt{t) : £ S X) is positive and 
bounded for all t. According to Eq. (j2.12p . this implies s^{t,to) ~ sinC(t,to) as |£| +oo, where 
C{t,tQ) is a nonzero function which depends on the system and whose form we do not need to 
specify. This is in conflict with the asymptotic behavior of s^(t,to) for the systems under study, 
given by s^(t, to) ~ as \i\ +oo for all (t, to) e I x I . 



C. Semiclassical states 



The explicit expression of the quantum unitary evolution for the single harmonic oscillator as a 
product of unitary operators (see t/ieoremO in sectioning turned out to be very useful to construct 
semiclassical states for some relevant one-dimensional dynamical systems. In particular, as stated 
in (I3.19p . the operator Tp{t) transforms the Lewis invariant ()3.18p into the time-independent free 
Hamiltonian (j3.17p . However, as we have shown in subsection \ VE\ there are obstructions that arise 
when dealing with systems of infinite oscillators -particularly, the possible nonunitarity of Tp{t)-, 
making the application of the techniques developed in section [7Vl particularlv difficult. In order 
to avoid these difficulties, we will probe an alternative procedure to construct semiclassical states 
that takes advantage of the unitary implementability of the quantum time evolution. We start by 
constructing the analogs of the minimal wave packets of the one-dimensional harmonic oscillator. 
Given a square summable sequence z = [z^ : ^ G X) E ^^(C), consider the state 

$(z) ._ g-|!z||2/2 ^-0 G , (5.19) 
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where the vacuum state corresponds in this context to z = 0, and ||z|| = Xl^ex \ ■ Vectors 
defined in this way appear as coherent superpositions of states with arbitrary number of particles. 
We then introduce the annihilation and creation operators in the Heisenberg picture corresponding 
to evolution backwards in time, 

ai{to,t) := U{t,to)aiU~^{t,to) = Ai{t,to)ai- Bi{t,to)a} , 
a}{to,t) := Uit,to)a}U-\t,to) = -Be{t,to)ai + Ai{t,to)a} , 

satisfying the Heisenberg algebra for all (t,to) & I x I ■ Here, A£(t,tQ) and B£{t,tQ) are the 
Bogoliubov coefficients defined in (|5.9p and (|5.10p . respectively. We then evolve the states ()5.19p 
in the Schrodinger picture, obtaining 

$(^)(t,to) :=C/(t,to)^^"^ = e-ll^ll'/2;7(to,t)exp(^z,a|)vI/o 

= e-il^ll V2 exp ( ^ z,a} (to, t)) (t, to) , 

with ai{to,t)<^'^'\t,to) = z^^>(^)(t,to), G X, and $(°)(t,to) = U{t,to)^o. By definition, the 
one-parameter family of states obtained in this way verifies the Schrodinger equation with initial 
condition and is closed under time evolution as well, U{t2,ti) ^^^\ti,to) = ^^^\t2,to) . The 
states ^^^\t,tQ) satisfy the properties stated in the following theorem. 

Theorem 6. Let z = {zi : £ G X) G £^(C) and to £ I- The position and momentum expectation 
values in the state <l>*^^)(t, to) = U{t,to)^^^^ satisfy 

qm{tM = («>(^)(t,to) |Qf $(^)(t,to)> = 2Q(t,to)Re(a^z^) + 2s<.(t,to)Re(^^z^), 
Pm(t,to) = {^'^'\t,to)\Pe^'^'\t,to)) = 2ce{t,to)Re{aeZe) + 2se{t,to)Re{peze), 

where {qHe,PHe) is the classical solution 113. S^) determined by the Cauchy data {qe,pe) = 
[2Re{a£Z£),2Ke{PiZ£)) at time to- Moreover, the mean square deviations of the position and mo- 
mentum operators with respect to the evolved state ^^^\t,tQ) satisfy 

^'S>(^){t,to)Q(- = \o^e^e{'t,to) + Q;^-B^(t,to)| , A^(z)(^^(p)P^ = \PiAi{t,to) + l3£B£{t,to)\ , 
irrespective of z. 

Remark 7. Let us consider the Gowdy models. As a consequence of the unitary implementability 
of the dynamics (jS.lip . we have B(t,to) G ^^(C) and, hence, |74^(t, to)| ~ 1 as \i\ — > +oo. We then 
obtain (see Fig. [3]) 

1 jW 
^$(-){t,to)<3^ ~ = ' ^'S>(-){t,to)^i^\(^A = \Hr when \l\ ^ +oo . (5.20) 

For fixed values of to, these asymptotic behaviors converge uniformly in t for time intervals away 
from the classical singularities. Note that these behaviors are the same that one would have 
expected to obtain if it had been possible to suitably extend the study developed in section [72 to 
field theories. We then conclude that the ^^^\t, to) vectors are coherent states far enough from the 
singularities. For the 3-torus model, the sequence z G ^^(C) is subject to satisfy a global constraint 
remaining on the system, given by X]^ez\{o} ^l-^^l^ = 0. In the Minkowskian case, expressions 
dOOH are valid for alU G Z \ {0} . 
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FIG. 3: Asymptotic behavior of the variances of the field and momentum operators Qi and for the states 
(t, to)) to < t & (0,7r), at high frequencies (hmit i +00) for the 3-handle and 3-sphere Gowdy models. 
Note the purely quantum behavior of these states as t approaches the classical singularities. These graphics 
can be considered as the limit of Fig. [5] for large frequencies. The asymptotic behavior for the 3-torus case 
is similar to the one represented in Fig. [2 identifying to = \i\ for each mode. 



VI. CONCLUSIONS 

In this paper, we have revised the unitary implementability of the quantum dynamics of a 
time-dependent harmonic oscillator (TDHO) and used the theory of invariants in order to define 
suitable semiclassical states for some relevant systems, such as the vertically driven pendulum 
or Gowdy-like oscillators. In particular, we have analyzed some important issues related to the 
associated Ermakov-Pinney equation, clarifying the need to introduce it as a natural way to obtain 
an evolution operator valid for all values of the time parameter. We must emphasize again that 
other elections different from the auxiliary Ermakov-Pinney equation may be problematic because 
of the singular behavior of the resulting evolution operator. We have shown that the Feynman 
propagator, usually derived by making use of more complicated path-integration techniques, can 
be obtained in a straightforward way within this scheme. The resulting formula has then been 
applied to calculate transition amplitudes, the instability of the vacuum state, and semiclassical 
states. Most of the calculations regarding the quantum evolution can be performed just by taking 
the classical dynamics into account, except for the presence of a phase that depends on the election 
of the quantum Hamiltonian. Nevertheless, this phase is irrelevant to answer all relevant physical 
questions such as the calculation of probability amplitudes or the evolution of quantum observables. 
It is important to remark that this phase, in contrast with the situation for the well-known TIHO 
system, cannot be eliminated in all cases by considering normal ordered Hamiltonians. 

Although we have concentrated our discussion on the quantum TDHO, our results can be easily 
extended to another interesting cases. For instance, our study is also applicable to the harmonic 
oscillator driven by an external, nonstationary, perturbative force characterized by a linear term 
in the position operator, 

H{t) = l{p^ + ^(t)Q') + fm , f^c\i). 

Indeed, this Hamiltonian can be transformed into the Hamiltonian ()3.9p just by introducing an 



/-dependent Glauber displacement operator 2l|, |58| 



Finally, we have extended the study of the unitary evolution of a single quantum harmonic 
oscillator to systems of infinite number of uncoupled oscillators with time-dependent frequencies, 
concretely, to the quantum field theory in Minkowskian space and the Gowdy cosmological models, 
providing a rigorous definition of the propagator. Here, the impossibility of unitarily implement- 
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ing some symplectic transformations turns out to be an obstacle to generalize the construction 
of semiclassical states through the eigenstates of the Lewis invariant [9|. Nevertheless, we have 
shown that the unitary implementability of the dynamics in appropriate Schrodinger representa- 
tions allows us to define suitable semiclassical states for these systems. In the case of the Gowdy 
cosmologies, they can be used to probe the existence of large quantum gravity effects in several 
ways. For instance, one may construct suitable regularized operators to represent the (three- or 
four-dimensional) metric of these models b y using arguments similar to those employed in the 
linearly polarized Einstein-Rosen waves 59.l6ol. 61 1 and the Schmidt model [l^. Calculating the 
expectation values of these operators in the coherent states, one may deduce the additional con- 
ditions (if any) that the sequences z S £^(C) (see theorem\^ in section\^ should satisfy in order 
to admit an approximate classical behavior. It is also important to analyze if the metric quantum 
fluctuations are relevant for all states. 

In addition, one may proceed as in [63] by appropriately promoting the quadratic invariant 
^"'^ Rabcd^"^^ R"''^'^'^ into a quantum mechanical operator. According to that reference, one should be 
able to unambiguously fix the operator order by requiring that the expectation values of this quan- 
tity in the coherent states exactly reproduce the classical results far from the singularities. In 
analogy with the results of [s^], even if the expectation values in other states (such as linear com- 
binations of coherent states) give nonclassical results, it is expected that the classical singularities 
persist in all cases. This physical consideration is supported by the purely quantum behavior of the 
uncertainties of the field and momentum operators in the coherent states at the classical spacetime 
singularities. 

The most natural way to extend the analysis developed in this article consists in considering 
generic nonautonomous quadratic Hamiltonians which contain the time-dependent harmonic os- 
cillators as particular cases. These systems can be analyzed from the perspective of some recent 
works on this subject (see 6^, 64, 6^]) in which Lie systems in quantum mechanics are studied from 
a geometrical point of view, developing methods to obtain the time evolution operators associated 
with time-dependent Schrodinger equations of Lie-type. These techniques may be successfully ap- 
plied to infinite-dimensional quadratic Hamiltonian systems by following a functional description 
similar to the one performed in this article. In particular, the different resulting factorizations for 
the time evolution operators may be especially useful to define alternative families of semiclassical 
states for these systems. 
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APPENDIX A: FACTORIZATION OF THE EVOLUTION OPERATOR 

In this appendix, we briefly summarize the construction of a well-defined factorization for the 
evolution operators of Minkowskian free fields and Gowdy cosmologies. The reader is referred to 
section [2 to revise the notation and definitions. A particularly useful way to proceed is to factorize 
U{t,to) as 



U {t, to) = Vp{t, to)np{t, to)Sp{t, to) , 
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with 



Vp{t,to) := Dp\t)Dpito), 
Sp{t,to) := Dp\to)Sp\t)Tp{to) , 
np{t,to) := T.\to)Rp{t,to)Tp{to) , 



where Dp[t) and Sp{t) are displacement and squeeze operators of the type defined in subsection 
\IIIE\ in such a way that 

V {t,to)aeVp{t,to) = [l + i\ai\ — TTT + ^TV ?7T ' (^^> 

V \Pdt) Pi{to)JJ \pe{t) pe{to)J 

V P^W pSo) Pe{io) \Pe{to) Pe{t) J J 

+ lai Oil—— - f3i]\ — — , A2) 

np^{t,to)ainp{t,to) = ^cos (^J^^ -^^^ +il^{ail3^ 

V Piiio)J J \Jto PlKV 

+ i \ 2aif3ipi{to)pe{to) - aj ( pj{to) + -oTTT ) ~ PePlito) 



Pi (to) 



* dr 



to 



X sin ( / ^ ) «^ (A3) 



and similarly for a| . Here, the solutions pi to the EP equations are conveniently selected as follows. 
For the Minkowskian free fields we choose pe{t) = l/y^j£\; for the 3-torus Gowdy model, we take 

Piit) = yy(Jo2(|^|t)+i^o'(W > (A4) 
whereas for the 3-handle and the 3-sphere Gowdy models we choose 



Piit) = fvr^Kcost) + ^^jicost) ] . (A5) 



Solutions (|A4p and (jASP have the asymptotic expansions 

peit) = + 0(|£r3/2) , p,(t) = C(t)l\tfl^ + 0(K|-^/2) , 

as \l\ — > +00 . Here, C(t) is a function of time which depends on the spatial topology and whose 
form we do not need to specify. It is then straightforward to check the unitary implementability of 
the transformations (|Aip , ()A2p , and (jASP in the Hilbert space by following the arguments employed 
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